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BASIC APPLICATION PROGRAMS 


The programs contained in this series are designed to demonstrate 
how the computer can be applied in a meaningful way to problems 

of many disciplines. The problems and the corresponding programs 
are, for the most part, quite simple and are designed to be 
"jumping off points" for students from the high school level on up. 


All of the programs, with very few exceptions, may be run on 
EduSystems 10 through 80, PDP-8 Family and PDP-11l Family computers. 
Exceptions are noted in the program descriptions. 


DEC welcomes contributions of programs and write-ups for use in 
this series. Contributors will be credited as program "Source". 
Please send programs and descriptions to Educational Marketing (5-2). 
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SCIENCE 


Demonstrates evolutionary mutations of 
pepper moths. 


Calculates and plots gas volumes at 
various pressures. 


Determines value of an unknown inductor. 
Tests knowledge of kinematics. 


Simulates the landing of a lunar module on 
the moon. 


Problems using Newton's second law. 
Demonstrates photo electric effect. 


Computes equivalent resistance for 
resistor networks. 


Produces tables of electric potential. 


AREA: Biology NUMBER : NAME: EVOLU 
SOURCE: Huntington Computer Project LANGUAGE: BASIC 


DESCRIPTION: The EVOLU program demonstrates the effects of genetic 
mutations of a sample of 10000 pepper moths. These mutations cause a 
small number of these moths to become somewhat darker in tone. However, 
it takes many years for this effect to show up, making it impossible to 
show in a school lab. In addition, the program allows the student to 
vary the mutation rate and to change the environment midway in the 
experiment. 


210 PRINT 
215 PRINT "ENTER RATE OF MUTATION (1 TO 140)" 
220 INPUT M 7 

250 LET Pa=100a0 

260 LET 7=Pa 

450 PRINT "YEAR IN WHICH ENVIRONMENTAL CHANGE OCCURS": 
460 INPUT X | 

464 LET L=1 

466 LEFT D=2 

470 PRINT "CHANGE FAVORS LIGHT MOTHS OR DARK CENTER 1 OR 2); 
480 INPUT E 

580 PRINT 

582 PRINT “YEAR"s "DARK MOTHS"»"LIGHT MOTHS" 

S524 PRINT 

610 FOR T=1 TO 3a 

620 IF T>=X THEN 65@ 

620 LET P1=0 

640 GO TO 710 

650 IF F<>2 THEN 630 

660 LET P1=INTCP1+e91*M*PQ+.5) 

670 LET PO=INT(Z7-P1+.5) : 

680 IF P1<Z% THEN 7)8 

69% LET P1=z 

700 LET Pasa 

T1A PRINT TsPlsPA 

720 NEXT T 

999 END 


READY 


USAGE: Type RUN. When asked, enter the mutation rate (from 1 to 10), 
the year in which the environment is to change (3 to 10) and whether 
the environmental change will favor light (type a 1) or dark ones 
(type a 2). The program then prints out the course of evolution over 


a thirty year period. V4 


FNTER RATE OF MUTATION (€1 TO 10)? 2 
YEAR IN WHICH ENVIRONMENTAL CHANGE OCCURS? 3 
CHANGE FAVORS LIGHT MOTHS OR DARK CENTER 1! OR 2)? 2 


YEAR DARK MOTHS LIGHT MOTHS 
1 A 10000 
2 4) | 1AgAA 
3 2A 98BAA 
4 396 9604 
Ss 5&8 9412 
6 776 9224 
7 964 9040 
8 1141 8859 
9 1318 8682 
10 1492 8SA6 
Le 1662 8338 
1? 1829 8171 
13 1992 BAAS 
14 2152 7848 
LD 23A9 7691 
16 P2463 7537 
17 2614 7386 
18 2762 7238 
19 L9IAT 71A93 
en 3A49 4951 
21 3188 6812 
22 3324 6676 
ne 3458 6542 
24 3589 6411 
as 3717 6283 
a6 3843 6157 
a7 3966 6034 
28 4A87 9913 
29 42065 9795 
34 4321 5679 


EDUCATIONAL BENEFITS: Where else can a student observe a thirty year 
reproductive cycle and watch the effects of mutations? More important, 
where else can the student repeat the cycle, trying new parameters and 
conditions? 


Many other common science experiments take a long time to run, or require 
the student to run it several times. To get meaningful results, they 
may often by simulated on the computer. 


AREA: Chemistry NUMBER : NAME: GASVOL 
SOURCE: NREL/SDC/DEC LANGUAGE: BASIC 


DESCRIPTION: The volume of a gas varies directly with the absolute 
temperature T (Kelvin) and inversely with the pressure P. If a certain 
quantity of gas occupies 500 cubic feet at a pressure of 53 pounds per 
square foot and an absolute temperature of 500 degrees, what volume will 
it occupy at 600 degrees absolute temperature and pressures from 100 to 
1000 pounds per square foot in increments of 50 pounds? 


The original conditions are used to solve for the constant K in line 
10. (K = VP/T). Then the new volumes are computed for varying pressures 
with T = 600° using the formula V = kKT/P. 


In the second part of the program lines 70 - 130 have been added to give 
a plot of the volume for the various pressures. 


How would you modify the program to deal with a more general case 
(1.e., other gasses and different temperatures) ? 


1A LET K=5AA*53/5AaM 

24 PRINT "A GASs AT 600 DEGREES KELVINs AND PRESSURE AS" 
25 PRINT “LISTEDs OCCUPIES THE VOLUME SPECIFIED." 

3A PRINT 

4A PRINT “PRESSURE, "VOLUME" 

5A FOR P=1AA TO 1AAH STEP 5A 

6A LET V=K*60A/P 

7A PRINT PsV 

8A NEXT P 

99 END 


C 


‘A GASs 


RUN 


A GASs AT 6A@ DEGREES KELVINs 


LISTED, 


PRESSURE 
1AA 
15a 
LAA 
25a 
SAA 
354 
AAD 
ASA 
SAA 
55a 
600A 
65a 
TAD 
TSA 
RAM 
RSA 
IAB 
954 
1AAA 


READY 


AND PRESSURE AS 


OCCUPIES THE VOLUME SPECIFIED. 


VOLUME 


7A PRINT P3 


RA LET X=INTCV*O/ 14465) 


318 

212 

159 
127-2 
106 
94-¢85714 
719¢5 
7A«66667 
A3-¢6 
57°¢81818 
a9 

48 «923088 
45242857 
AP 
39275 
37¢41176 
3533333 
33047368 
3168 


90 FOR L=1 TO xX 
1A4 PRINT " "5 


110 NEX 


TL 


120 PRINT "x" 


130 NEX 
999 END 


READY 


RUN 


LISTED, 


PRESSURE 
100 
150 
200 
2590 
300 
350 
AGA 
450 
5A 
559 
6AM 
650 
700 
750 
BAO 
R50 
900 
950 
1000 


T P 


AT 600 DEGREES KELVINs 


AND PRESSURE AS 


OCCUPIES THE VOLUME SPECIFIED. 


PLOT OF VOLUME 


re me 


AREA: Electronics NUMBER: NAME: INDUC 
SOURCE: Digital Equipment Corp. LANGUAGE: BASIC 


DESCRIPTION: This program combined with the indicated laboratory 
measurements will determine the value of an unknown 
inductor, such as a choke or transformer winding. 


APPROACH: The impedance of an inductor 1s made up of two 
quantities; inductive reactance (X_) and d.c. re- 
Sistance (R). First, measure the vesistance with 
an ohmmeter. 


Next, determine the voltage (E) and current (I) of 
an a.ec. circuit where the unknown inductance con- 
stitutes the entire load. Ohm's Law may be applied 
to find the impedance, Z = B/I. 


We can now separate R from X_ by means of a vector 
diagram in ge R becomes the base and Z, the hypo- 
tenuse. Once has been obtained, the inductance L 
may be found trl the formula: 


In the formula, f is the frequency. A simple com- 
puter program can replace the various calculations 
and vector diagram. The program plus a sample run 
are shown on the back. 


Ah 
NMILLIAMME TER 


VARIAC 4 UNKNowN 
oR LINE INDUCTOR - 
VOLTAGE Z= 2000 OHMS 





XL = 1830 OHMS CAPPROK) 


5 PRINT “PROGRAM INDUC" 
6 PRINT 
1A PRINT "TO FIND THE VALUE OF AN UNKNOWN INDUCTOR," 
11 PRINT “INPUT THE FOLLOWING QUANTITIES:" 
12 PRINT 
15 PRINT “RESISTANCE"S 
16 INPUT R 
PA PRINT "“VOLTAGE"S 
21 INPUT FE 
25 PRINT "CURRENT" 
26 INPUT I 
3A PRINT "FREQUENCY"''S 
31 INPUT @ 
35 LET Z=E/I 
4A LET XK=SQRCZ*Z7-R*¥R) 
45 LET L=X/7C2*3-14159*0) 


5A PRINT 
51 PRINT “INDUCTANCE IS"L"HENRIES-" 
99 END ‘ 

READY 

RUN 


PROGRAM INDUC 


TO FIND THE VALUE OF AN UNKNOWN INDUCTOR; 
INPUT THE FOLLOWING QUANTITIES: 


RESISTANCE? 880A 
VOLTAGE? 1AaA 
CURRENT? «A5 
FREQUENCY? 6A 


INDUCTANCF IS 4862268 HENRIES.- 


READY 
RUN 
PROGRAM INDUC 


TO FIND THE VALUE OF AN UNKNOWN INDUCTORs, 
INPUT THE FOLLOWING QUANTITIES: 


RESISTANCE? 4a4@ 
VOLTAGE? 115 
CURRENT? e«@46 
FREQUENCY? 64 


INDUCTANCFE IS 6¢546A29 HENRIES- 


READY 


AREA: Physics NUMBER : ; NAME: KINEMA 
SOURCE: Digital Equipment Corp. LANGUAGE: BASIC 


DESCRIPTION: The KINEMA program tests the student's knowledge of kine- 
matics. A ball, or any other object, is thrown up at a certain velocity 
(in meters per second.) The student is then asked to calculate certain 
information about the resulting trajectory. The computer then checks the 
response to see if it is within 15% of the correct answer. 


180 PRIN? 

105 PRINT 

186 LET 0=0 

110 LET V=5+INTC35*RNDCO)) 

111 PRINT "A BALL IS THROWN UPWARDS AT"V"METERS PER SECOND" 
112 PRINT 

115 LET A=e@5*Vte 

116 PRINT "HOW HIGH WILL IT GO CIN METERS)"; 
117 GOSUB 560 

129 LET A=V/5 

122 PRINT “HOW LONG UNTIL IT RETURNS CIN SECONDS )"'s 
124 GOSUB 5A0 

130 LET T=1+INTC2*V*RNDC1)2/10 

132 LET A=V-10*T 

134 PRINT “WHAT WILL ITS VELOCITY BE AFTER"T*'SECONDS"s 
136 GOSUB 506 

140A PRINT 

150 PRINT Q"RIGHT OUT OF 3-6"; 

160 IF Q<2 THEN 100 

170 PRINT * NOT BAD." 

18@ GO TO 146 | 

500 INPUT G | 

502 IF ABSC(CG-A)/A)<-15 THEN 510 

504 PRINT “NOT EVEN CLOSEeeeoe*’ 

586 GO TO 512 

510 PRINT “CLOSE ENOUGH." 

511 LET Q=0Q+1 

512 PRINT “CORRECT ANSWER ISA 

520 PRINT 

530 RETURN 

999 END 


READY 


USAGE: Load KINEMA, then type RUN. Then try to answer the questions. 
Distances should be entered in meters, time in seconds, and velocity 
in meters per second (plus for up, minus for down.) 


RUN 


A BALL IS THROWN UPWARDS AT 13 METERS PER SECOND 


HOW HIGH WILL IT GO CIN METERS)? 8 
CLOSE ENOUGH. 
CORRECT ANSWER IS 845 


HOW LONG UNTIL IT RETURNS CIN SECONDS)? 3 
NOT EVEN CLOSEeccee 
CORRECT ANSWER IS 26 


WHAT WILL ITS VELOCITY BE AFTER 1+¢7 SECONDS? -4 
CLOSE ENOUGH. 
CORRECT ANSWER [5-4 


® RIGHT OUT OF 3- NOT BAD- 


A BALL IS THROWN UPWARDS AT 30 METERS PER SECOND 


HOW HIGH WILL IT GO CIN METERS)? 590 
CLOSE ENOUGH. 
CORRECT ANSWER IS 45 


HOW LONG UNTIL IT RETURNS CIN SECONDS)? 6 
CLOSE ENOUGH. 
CORRECT ANSWER IS 6 


WHAT WILL ITS VELOCITY BE AFTER 2«¢8 SECONDS? 4 
NOT EVEN CLOSEceee 
CORRECT ANSWER IS 2 


2 RIGHT OUT OF 3- NOT BADe 


EDUCATIONAL BENEFIT: The key educational benefit here is the powerful 
motivating force that the computer has. The work that the student is 
being asked to do is really no different than rote homework exercises. 
But the computer makes it a challenge. It is very frustrating to have 
a dumb computer keep telling you you're wrong. Try it yourself and see 


if it doesn't make you want to brush up on your physics. That's real 
motivation. 


Fe Dw 3 
AREA: Games/Physics NUMBER: NAME : LUNAR 
SOURCE: Digital Equipment Corp. LANGUAGE: BASIC 


DESCRIPTION: The program represents an exact simulation of an Apollo lunar 
landing module during the final descent. This portion of the descent would 
normally be controlled by the on-board computer backed up by a computer 
located on Earth. However, to exercise your knowledge of physics (and 

make an interesting game), both computers simultaneously have had a 
malfunction; therefore, you are on your own to safely land the spacecraft. 


To make a soft landing, you may reset the burn rate of the retro rockets 
every ten seconds. You have a choice of not firing at all (burn rate=0) 
or firing at a rate between 8 and 200 lbs. per second. You have 16500 
lbs. of fuel. If the rockets were not fired, your estimated free fall 
impact time is 120 seconds. The capsule weight is 33,000 pounds 


USAGE: This program runs on EduSystem 20 and larger systems. Type SCR 

to get rid of any existing programs and load LUNAR. Type RUN to commence 
landing procedure. When a "?" is typed under burn rate, enter your burn 
rate and hit the return key. Remember, the only acceptable burn rates are 
O or any number between 8 and 200. 


After you're down and the computer has typed "READY", if you wish to try 
it again, simply type RUN. Good luck! 


RIN 

GROUND CONTROL CALLING LUNAR MODULBececece 

ON-ROARD AND GROUND COMPUTERS HAVE SIMULTANEOUSLY 
MALEF'TINCTIONFD CTHEY WEREN'T DEC MACHINES) 

MANUAL CONTROL IS NECESSARYe 


CAPSULE WEIGHT 325590 LBS AVAILABLE FUFL 1635990 LBS : 
FSTIMATED FREE FALL IMPACT TIME 1298 SECONDS 


SET RETRO ROCKET BURN RATE EVERY 10 SECONDS TO ANY VALUE ° 
BETWEEN @ LBS/SFC CFREE FALL) AND 
288 LBS/SEC CSTRONG BRAKING) 


GOOD LUCK ri? 


SEC MI + FT MPH LB FUEL BURN RATE 

a i129 @ 3699 16590 2 o 

19 199 5916 3636 16500 2 0 

20 99 4224 3672 16509 2 0 

30 89 2994 3798 16500 2 0 

Ag 79 1956 3744 16500 2 

50 68 3969 3780 16580 2 Q = 
69 58 1056 3816 16500 20 | 
70 AT 2904 3852 16500 2 200 — 
89 37 1884 3482-868 14500 2 209 

99 28 1191 3086+ 708 12500 2 200 

189 29 1251 26592654 19500 ? 209 

110 13 2549 21960947 8500 2 200 

120 8 370 1692-634 6500 2 200 

130 4 658 11392138 A598 2 200 

149 1 4294 526-5976 2500 2 109 

159 Q A942 212.240 1590 2 45 

160 ®@ 1864 84018261 1950 2 20 

17D 9 909 45-91246 850 2.17 

180 B 438 18.10655 680 2.13 

199 B 269 42886334 550 2 19 

250 ® 213 20768577 450 2 10 : 
210 % 189 -4252781 350 2 902 

220 @ 179 29479145 258 2 903 

230 ® 165 » 8864265 165 2 904 ; 
240 @ 157 © 2329987 71 2 Tel ; 


RAN OUT OF FUEL AT 258 SEC 

ON THE MOON AT 254-9616 SEC - IMPACT VELOCITY 23-8930908 MPH 
CRAFT DAMAGEe YOU'RE STRANDED HERE UNTIL A RESCUE MISSION 
ARRIVES. HOPE YOUR OXYGEN HOLDS OUT! 


REA DY 


READY 


S 
Uv 
7 


ORYAUNEDWO- 


«192 
45183 
200119 
275311 
wWr12 


PRI "GROUND CONTROL CALLING LUNAR MODULFeeeece"™ 

PRI “ON-BOARD AND GROUND COMPUTERS HAVE SIMULTANEOUSLY" 

PRI "MALFUNCTIONED CTHEY WEREN*T DEC MACHINES)" 

PRI "MANUAL CONTROL IS NECESSARY." 

PRI 

PRI “CAPSULE WEIGHT 325,500 LBS AVAILABLE FUEL 1625998 LBS" 
PRI “ESTIMATED FREE FALL IMPACT TIME 128 SECONDS" 

PRI 

PRI “SET RETRO ROCKET BURN RATE EVERY 10 SECONDS TO ANY VALUE” 
PRI “BETWEEN @ LBS/SEC CFREE FALL) AND" 

PRI "20990 LBS/SEC CSTRONG BRAKING)" 

PRI 

PRI "GOOD LUCK rie 

PRI 

PRI “SEC's*MI + FFTs "MPH">"LB FUEL's"BURN RATE" 

PRI 

A=120 

V=1 


> M=33909 


N=165390 

G=e4G1 

Z=128 

PRI Ls INTCAIs INTC 5289*CA-IN TCAD) s 369G0*Us MaNs 
INPUT kK 

T=193 

IF M-N<-e991 THEN 72 

IF T<e991 THEN 3903 


S=T 

IF M>=N+S*K THEN 50 
S=(M-N)/K 

GOSURB 173 


IF I<=8 THEN 139 

IF V<=8 THEN 69 

IF J<9 THEN 159 

GOSUB 1293 

GOTO 49 

PRINT "RAN OUT OF FUEL AT'"L"SEC"™ 
S=(€-U+tSARCU*XU+F2*A*G)I/G 

V=V+GeS 

L=L+S 

W=3690*V 

PRINT "ON THE MOON AT'L"™SEC - IMPACT VELOCITY" W"MPH" 

IF W>le2 THEN 92 

PRI “PERFECT LANDING! CLUCKY)" 

STOP 

IF W>1@ THEN 19@ 

PRI "GOOD LANDING CCOULD BE BETTER)" 

STOP 

IF W>69 THEN 119 

PRI "CRAFT DAMAGEe YOU'RE STRANDED HERE UNTIL A RESCUE MISSION" 
PRI "ARRIVESe HOPE YOUR OXYGEN HOLDS OUT!® 

STOP 

PRI “SORRYs BUT THERE WERE NO SURVIVORS -- YOU BLEW ITe"™ 
PRI “IN FACTs YOU BLASTED A NEW LUNAR CRATER'W*-e2777"FT DEEP" 
STOP 


2 129 
2vw 121 
228 128 
239 1283 
255 1284 
240 185 
248 136 
79 131 
255 132 
285 141 
270 142 
215 158 
290 151 
285 152 
7q0 («168 
2s (161 
xm 162 
PS 163 
3° 164 
as 172 
e796 «171 
ez 188 
3% 181 
2A 


REANY 


qpisut 
oun wa 
1 + 
NNN 

% 

x 


ay eS 


RET 

IF S<e905 THEN 89 

D=V+SQRCV4V+t2#A*C G-Z*K/M) ) 

S=2%xA/D 

GOS 170 

GOS 129 

GOTO 139 

W=C1-M*G/CZ*K))/2 

S=M*U/ CZ*K* CW+ SOR CWHWtUSZI ID +285 

GOS 172 

IF I1<=9 THEN 139 

GOS 12@ 

IF J>@ THEN 49 

IF V>@ THEN 150 

GOTO 48 

Q=S*K/M 
J=UtEKS-Z*OKC1L4+0% Ce 540401 734+0%0C2254+0/5)))) 
T=A-GRS*SSO-URS+ZESHKOR Ce 54+0*C1/64+0*01712840/20))) 
RET 

END 


AREA: Games/Physics NUMBER: NAME: LUNAR2 
SOURCE: Digital Equipment Corp. LANGUAGE: BASIC 
(EduSystem 10 only) 


DESCRIPTION: This program is a compressed version of the lunar module 
landing program originally written in 8K FOCAL. The program is also known 
as LUNAR, ROCKES, and ROCKET. 


The program represents an exact simulation of an Apollo lunar landing 
module during the final descent. This portion of the descent would 
normally be controlled by the on-board computer backed up by a computer 
located on Earth. However, to exercise your knowledge of physics (and 
make an interesting game), both computers simultaneously have had a 
malfunction; therefore, you are on your own to safely land the spacecraft. 


To make a soft landing, you may reset the burn rate of the retro rockets 
every ten seconds. You have a choice of not firing at all (burn rate=0) 
or firing at a rate between 8 and 200 lbs. per second. You have 16500 
lbs. of fuel. If the rockets were not fired, your estimated free fall 
impact time is 120 seconds. The capsule weight is 33,000 pounds. 


USAGE: This program runs on EduSystem 10 (4K) BASIC only. Retain only 
the SQR and RND functions (Option F). Type SCR to get rid of any exist- 
ing programs and load LUNAR2. Type RUN to commence landing procedure. 
When a "?" is typed under burn rate, enter your burn rate and hit the 
return key. Remember, the only acceptable burn rates are O or any 
number between 8 and 200. 


After you're down and the computer has typed "READY", if you wish to try 
it again, simply type RUN. Good luck! 


LIST 


PRINT “SEC’s"MI + FT's"MPH's"LB FUEL's "BURN RATE’NPRINT 
A=120NV=1NM=33000NN=1 6500NG=1 -OAQ0000E-O3NZ=1-8 

PRINT Lo INTCA)3 INTCS5280* CA-INTCA)) 5 3609* VsM-Nes\NINPUT KNT=1QNPRINT 
IF M-N<-@01 THEN 4INIF T<-001 THEN 21NS=TNIF M<N#S*K THEN S=(M-NO/7K VO, 
GOSUB 91 NIF I<=@ THEN 71 NIF V<=@ THEN 38 NIF J<@ THEN 81 | 
GOSUB 61 NGO TO 31 

PRINT “FUEL OUT AT'SL3"SEC"NS=C- V+ SQRCV* V+ 2*A* G) 07 GN V= Vt G¥ SNL=ELtS 
W=3600*VNPRINT “ON MOON AT’SLS""SEC - IMPACT VELOCITY'S Ws "*MPH™ 

IF W>1@ THEN 56 NPRINT “GOOD LANDING"\ STOP 

IF W>6@ THEN 58 NPRINT “BAD JOLT’"NSTOP 

PRINT “NO SURVIVORS-BLASTED NEW LUNAR CRATER" W*-27773"FT DEEP"\STOP 
L=L+SNT=T- S\M=M- S¥KNA=INV=JNRETURN 

IF S<5-8@00000E-03 THEN 51 \D=V+SQRCV#Vt2*A* (C G=Z*K /M) INS=2*A/D 

GOSUB 91 NGOSUB 61 NGO TO 71 

W=C1-M* G/CZ¥*K I) /OENS=M*V/CZ¥EK* C Wt SQRC WH WtV/Z99094+-8S\NGOSUB 91 7 
IF I<=@ THEN 71 \NGOSUB 61 NIF J>@ THEN 31 NIF V>@ THEN 81 \GO TO 31 

Q= S*K /MNJ= V+ G* S+Z* C= Q- Q*0/72-Qt 3/3-Qt 4/ 4-Qt 5/5) 

T=A~ G¥S¥*S/2-V¥S+Z¥*S¥*CQ/2+ Qt 2/ 64+ Qt 3/124 Qt 4/26+ Qt 5/30) \ RETURN 

END 


aac MI + FT MPH LB FUEL BURN RATE 

B 128 @ 3600 16500 20 

10 109 5016 3636 16500 20 —_, 
20 99 seen 3672 1658 20 

30 $9 2904 3708 16580 208 

4B 79 1056 3744 16502 20 

58 68 3960 3780 16580 20 

60 58 1056 3816 1650G (2200 

70 48 S64 3446-868 14580 2200 

82 39 399 3058-708 12580 2200 

90 31 987 2623-654 12500 20 
190 23 4483 26590654 10500 2200 

110 17 50! 2196947 8500 20 

120 19 4975 2232-947 8500 2200 

130 5 2268 1728-634 6500 2200 a 
1 40 1 2028 1175-137 4500 2200 a, 
ON MOON AT 144-8262 SEC - IMPACT VELOCITY 887-6962 MPH 


NO 


SURVIVORS-BLASTED NEW LUNAR CRATER 246-4883 FT DEEP 


READY -e 


AREA: Physics NUMBER : NAME: NEWTON 
SOURCE: Project SOLO, Pittsburgh, Pa. LANGUAGE: BASIC 


DESCRIPTION: (This problem assumes the student has basic know- 
ledge of Newton's three laws) 


The first problem we will solve is as follows: 
A car at rest at a traffic light weighs 5.63 * 10 °newtons. 
(a) What is the mass of the car? 
When the light changes to green, it accelerates at 3.18 m/sec* 
(b) If friction is negligible, what must be the force between 
the rear wheels and the road? 
(c) What is the velocity of the car when it has moved 30 meters? 


Remembering that gravitational force is 9.8 m/sec*and W=MA, we 
can calculate the mass of the car directly in statement 20. We 
can then make use of the relationships: 


F = MA 
V Vo t+ A* Tor =Vo +2*A*D 


5S PRINT “INPUT Ws As D3 
10 INPUT WsAsD 
OA LET M=V/9eR 

3A LET F=Mx*xA 
AA LET V=2*AxD 
SA PRINT "M="M,"Pa"F, "yay 
99 END 


PFADY 
FUN 


INPUT Ws As D? Se63E35 3-185 3A 
M=z 574-4898 F= 1826-e878 V= 190-8 


If some other force in addition to the gravitational force acts 
on a falling body, it's acceleration will not be 9.8 m/sec ~. 
Consider the following problem: The Saturn V rocket has a mass 
of 3.18 x 10°kg. The thrust of its first stage engines is quoted 
as being 7.5 * 10° lbs. What is the initial acceleration of the 
rocket as it blasts off for the moon. 


In this case, Newton's second law is expressed: 
(F + W)= MA 


Since the direction of rocket acceleration is upward, W has a 
negative value. 


5 PRINT “INPUT Ms F's 
14 INPUT MsF 

PA LET W=-M*9.R 

3A LEFT A=(F+W)/M 

“4A PRINT "ACCEL ="'A 
99 FND 


RUN 


INPUT Ms F? 3e1R8EF6s 3¢34E7 
ACCEL = e¢7TA31447 


FOR MORE INFORMATION: The teaching module "Newton's Laws of Motion" 
containing six additional excellent problems is available complete 
from: 

Dr. Thomas Dwyer 

Project SOLO, Computation Center 

University of Pittsburgh 

Pittsburgh, Pa. 15213 


AREA: Electricity NUMBER: 8-601 NAME: Resist 
SOURCE: Computer Methods in Mathematics LANGUAGE: BASIC 


J DESCRIPTION: The problem is to computer the equivalent resistance 
between points A and B for the circuits below: 


RI 


R3 


<3 


APPROACH: We must remember that in series, resistance is additive, 
that iss: 


Rp = Ry + RQ --- 


— Whereas in parallel, the resultant resistance is com- 
puted as follows: 


Given these relationships, it is a simple matter to 
write a program to compute the resultant resistance for 
the two circuits above. 


5 PRINT “PROGRAM RESIST" 
6 PRINT 
16 PRINT "RI"; 
15 INPUT RI}! 
2A PRINT "RO" S 
25 INPUT R2 
Z 34 PRINT "R3"'$ 
35 INPUT R3 
4A LET R42#1/C1/7R1+1/7R2)+R3 
52 PRINT 
60 PRINT "R4 =""3R4 
99 END 


RUN 
PROGRAM RESIST 
R1? 330 


RE? 22d 
R3? 110A 


RUN 
PROGRAM RESIST 
R1? 10 


R2? 10 
R3? 1A 


READY 


For the second circuit, we need change only one statement 


in our program: 


40 LET R4=1/(€1/R1+1/R2+1/7R3) 


RUN 

PROGRAM RESIST 
R1? 330 

R2? 220 

R3? 110 


R4 = 60 


READY 

RUN 

PROGRAM RESIST 
R1? 140 

R2? 10 

R3? 10 


R4 = 3¢333333 


READY 


EDUCATIONAL BENEFITS: 


The computer can be used as a tool to aid in 
solving all kinds of practical problems. It 
is not necessary for a student to get bogged 
down in mathematical calculations when he is. 
eagerly learning about other things. 


PUN 
PROGRAM UFLEC 


WHAT IS THE VALUF OF 9€1)? 14 
WHAT IS THE VALUF OF Q¢€2)? 1! 


WHAT APF THE DISTANCE LIMITSs RC1)s RC2)s BETWEEN 
VHICH YOU DESIRE THE VALUFS TO RUN? 122 
HOW MANY SURDIVISIONS OF THE INTERVAL DO YOU WISH? 12 


DISTANCE» R ELFCTRIC POTENTIALs U 


METERS JOULES 

1 BeDBRGAGE +10 
1°¢AR3333 Be 296615E+10 
1-166667 Te TALAAGBE+10 
1225 7° 19A4O00E+194 
16333333 6° TAIAAAE+14 
12416667 6¢344471E+10 
1-5 5 -992DAAE+140 
1258232333 5 -6764A32E4+14 
1-666667 5 - 3928 AAE+140 
1-75 Se 1LSEAANE+14 
128323333 4e¢9A25A6E4+10 
12916667 4e6RIO3B91EF+1A 
? Ae AQD4AAAEF IO 

PF ANDY 

PTTIN 


PROGPAM UELEC 


VHAT 1S THF VALUF OF Q¢€1)? 14a 
WHAT IS THF VALUE OF @(2)? 5 


WHAT ARE THE DISTANCE LIMITSs RC€1)5 RC2)s BETWEEN 
VHICH YOU DESIRE THE VALUES TO RUN? e- 

MORE? e535 

HOV MANY SUBDIVISIONS OF THE INTERVAL DO YOU WISH? 9 


DISTANCEs R ELECTRIC POTENTIALs U 
METERS JOULES 


63) 


8 -9BBAAGE+11 
4-494A00E+11 
2 e996AG0GF +11 
2e2A4TAAPE+1 1 
1°79 7600EF+11 
1¢498AAAE+11 
1-28 4000E+11 
1-1235AK0E+11 
9-9R6667E+14 
8 e9BBAAAE+1IA 


® 
io 
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nn 
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READY 


AREA: Physics/Electricity NUMBER: NAME: UELEC 
SOURCE: J. Martin, Rutland High School LANGUAGE: BASIC 
DESCRIPTION: This iS a program which produces a table of values for 
tne electric potential energy, in joules, at various 
Separations, in meters (R), of two electric charges, 
Ql and Q2. 
Uist ac (K*Q1*Q2)/R 
The user inputs the values of Ql, Q2, Rl, R2 and the 
number of subdivisions desired in the total interval 
of distance R2-Rl. Note that Rl must be less than R2. 
5S PRINT "PROGRAM UFRLEC" 
6 PRINT 
100% PRINT “WHAT IS THE VALUF OF @(€1)"S 
110 INPUT Q@1 
120 PRINT “WHAT IS THE VALUE OF Q(€2)"'S 
130A INPUT @2 
140 PRINT 
150A LET K=8&e9R& EQ : 
1640 PRINT “WHAT ARF THE DISTANCE LIMITSs RC€1)5 RC2)s BETWEEN” 
170A PRINT “WHICH YOU DESIRE THE VALUES TO RUN"; 
184A INPUT R1»sR2 
190A PRINT "HOW MANY SUBDIVISIONS OF THE INTERVAL DO YOU WISH"’S 
CAA INPUT S 
2AS PRINT 
21A PRINT “DISTANCE, R"s"*ELECTRIC POTENTIAL, U" 
215 PRINT °° METERS*'s °° JOULES" 
216 PRINT 
2°2A FOR R=R1 TO RO STEP CR@-R1)/7S 
225 LET U=K*01*Q2/7R 
230A PRINT ° ak a "su 
25a NEXT R 
999 END 


READY 


AREA: Physics NUMBER: NAME: PHOTOE 
SOURCE: Huntington Computer Project LANGUAGE: BASIC 


DESCRIPTION: When light of short wave-length falls on a metal surface, 
electrons are ejected from the metal. According to the description of 
this phenomenon by Einstein, there is a maximum wavelength for each 
metal above which no electrons are emitted. This program lets the user 
find this critical wave-length. 


APPROACH: The program simulates the actual experiment, where the metal 


is placed in a vacuum and bombarded with soft X-rays. The number of 


electrons ejected is collected and counted by means of an ammeter. 

The program types out three ammeter readings for each wave-length. After 
this has been done, the intensity may be increased and the experiment 
rerun. The typed out values coule then be graphed to show the critical 
wave-length. 


To pick a metal enter its coefficient as follows: 


Silver - 308 
Bismuth -338 
Cadmium -318 
Lead - 340 
Platinum .385 


EDUCATIONAL BENEFITS: Students very rarely get to work with the same 
concepts that Einstein worked with, or indeed with any of the great 20th 
century work in physics. The reason is simple: modern physics uses 
laboratory gear that is way beyond the means of most schools. The con- 
cepts themselves are often easy enough for the student; it is only the 
lack of equipment that holds him down. As this program demonstrates, 
the computer can fill this gap by simulating equipment. And unlike most 
lab equipment which is used with one unit, then lies idle the rest of 
the term, a computer does all kinds of things all year long. Its 
flexibility means it is never idle. 


146 PRINT | | it ged ee 
110 PRINT “ENTER CO-EFFICIENT FOR METAL" 
120 INPUT VA 
460 LET K=INTC1+2*RNDCX)) 
470 PRINT 
480 PRINT * ">"*MEASURED CURRENT (MICRO-AMPERES)" 
ASA PRINT “WAVELENGHTH"s"TRIAL 1°%s"TRIAL 2"s"TRIAL 3" 
SAM FOR L=e42 TO -25 STEP -.A2 
S10 LET M=INTC1AGG/L) 
520A PRINT Ms | 
530A FOR J=1 TO 3 
540 IF L>V@ THEN 570 
550 LET I=SQORCINTC25*RNDCX))) 
S560 GO TO 580A. 
570 LET I=SQRCK*K*100+INTC35*RNDC(X) )?) 
S5SB8A LET N=INTC10*1+-5)/10 
590 PRINT No 
600 NEXT J 
610 PRINT 
620 NEXT L 
625 PRINT 
630 PRINT “INCREASE LIGHT INTENSITY (€1 FOR YES» 2 FOR NO)"'$ 
660 INPUT 92 
670 IF Q2=2 THEN 999 
690 PRINT "BY WHAT FACTOR (1 TO 1@)"3 
700 INPUT F 
720 LET K=K*F 
730 GO TO 47a 
999 END 


READY 


RUN: 


FNTER CO-EFFICIENT FOR METAL? +340 


MEASURED CURRENT (MICRO-AMPERES) 


WAVELENGHTH TRIAL 1 TRIAL 2 TRIAL 3 
2380 1965 11-2 10-5 
2500 145 10 10-8 
2631 108 19-8 —«*18 
2777 144 10-5 193 
2941 — 4e7 2.4 49 
3125 2.4 Q 46 
3333 48 4e5 2.6 
3571 363 3.3 2.2 
3846 Ad 4 AeA 


INCREASE LIGHT INTENSITY (1 FOR YESs 2 FOR NO)? 1 
BY WHAT FACTOR (1 TO 10)? S 


MEASURED CURRENT (MICRO-AMPERES ) 


WAVELENGHTH TRIAL 1 TRIAL 2 TRIAL 3 
2380 50 | 5Bo1 50-2 
2500 50-2 | 5G-2 50-2 
2631 50-2 5A. 58-3 
29777 50-2 58-2 50-3 
(2941 4e8 | 307 104 
3125 3-2 | 367 307 
ree AS | 2 : ; 124 4-6 
3571 ate Le? Ae fo BR —Me2 
3846 ee S. < « eS ‘ee Aol 


INCREASE LIGHT tNTENSITY (1 FOR YES» 2 FOR NO)? 2 
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Digital Equipment AB 
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Cable: Digital Stockholm 
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Digital Equipment Corp. A/S 
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Telephone: 02/68 34 40 
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Digital Equipment Aktiebolag 
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Digital Equipment Corporation S.A. 
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81 Route de |'Aire 
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Telephone: 42 79 50 Telex: 22 683 
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Scheuchzerstrasse 21 
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Digital Equipment Corporation Ltd. 
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Digital Equipment of Canada, Ltd. 
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Telephone: (613)-257-2615 TWX: 610-561-1651 


Telex: 25297 


Telex: 19079 DEC N 
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230 Lakeshore Road East, Port Credit, Ontario 
Telephone: (416)-274-1241 TWX: 610-492-4306 
MONTREAL 
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Digital Equipment Australia Pty. Ltd. 
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P.O. Box 491, Crows Nest 

N.S.W. Australia 3065 
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Ex. 2729 


Telex: AA20740 


Telephone: 63-1339 Telex: AA82825 
NEW ZEALAND 

Digital Equipment Corporation Ltd. 
AUCKLAND 


Hilton House, 430 Queen Street, Box 2471 
Auckland, New Zealand 
Telephone: 75533 


JAPAN 

Digital Equipment Corporation International 
TOKYO 

Kowa Building No. 17, Second Floor 

2-7 Nishi-Azabu 1-Chome 

Minato-Ku, Tokyo, Japan 


Telephone: 404-5894/6 Telex: TK-6428 


tts, U.S.A. 01754 * Telephone: From Metropolitan Boston: 646-8600 « Elsewhere: (617)-897-5111 
TWX: 710-347-0212 Cable: DIGITAL MAYN Telex: 94-8457 
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Virrey del Pino 4071, Buenos Aires 


Telephone: 52-3185 Telex: 012-2284 
BRASIL 

RIO DE JANEIRO — GB. 

Ambriex S.A. 
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CHILE 


SANTIAGO 
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INDIA 


BOMBAY 
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Mexico 12, D.F. 
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